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Abstract 

o ■ 

, The standard cold dark matter model with a cosmological constant is the best fit to cosmological 

observations and to galaxy rotation curves. However, unless there is a direct detection of dark 
matter, either in the laboratory or in astronomical observations, one should allow for modified 
gravity theories such as MOND or Scalar- Tensor- Vector Gravity [STVG] as possible explanations 
' for flatness of galaxy rotation curves. The STVG theory due to Moffat and collaborators modifies 

general relativity by the addition of a massive vector field, and the vector field coupling constant, 
its mass, and the gravitational constant, are dynamical scalar fields. The theory is shown to yield 
a modified acceleration law which has a repulsive Yukawa component added to the Newtonian law 
of gravitational acceleration, and which can explain the observed flatness of a large class of galaxy 
■ rotation curves by fixing the values of two free parameters [a mass scale and a length scale]. In 

| the present paper we provide a possible insight into the success of the STVG theory, by staying 

within general relativity, and by considering the effect of polarization [quadrupole approximation] 
on the averaged gravitational field inside a galaxy, due to the pull of neighboring galaxies. This 
effect is analogous to the polarization induced modification of averaged electromagnetic fields in an 
electrically charged macroscopic medium, and was studied by Szekeres in the context of propagation 
j—j ' of gravitational waves. The study was generalized to the case of a static weak-field approximation by 

C/3 , Zalaletdinov and collaborators, who showed that the effect of quadrupole polarization is to modify 

Poisson's equation for the gravitational potential to a fourth order [biharmonic] equation. We show 
that, remarkably, the biharmonic equation implies a modification of Newtonian acceleration which 
is precisely of the same repulsive Yukawa form as in the STVG theory, and the corrections could in 
^ ' principle be large enough to explain flatness of the rotation curves. 

in 
t> 

m ■ 1 Introduction 

Dark matter is the most well accepted explanation for the flatness of galaxy rotation curves, and it also fits 
very well into the standard ACDM paradigm for cosmology. However, unless there is a direct detection of 
one or more popular dark matter candidates in the laboratory, or via astronomical observations, alternate 
possible explanations should not be excluded. These alternatives include modified theories of gravity 
such as MOND (Modified Newtonian Dynamics) and the Scalar- Tensor- Vector Gravity theory developed 
by Moffat and collaborators [T], [2], [3]. It is the latter which we will be concerned with, in the present 
paper. STVG is a generalization of general relativity which includes a massive vector field, and also 
three dynamical scalar fields G{x), fx(x) and ui(x), which are respectively the gravitational 'constant', 
the mass of the vector field, and its coupling to gravity. It has been shown that STVG possesses a 
modified law of acceleration which can explain flat rotation curves for a large class of galaxies after fixing 
two parameters, a mass scale and a length scale, to universal values. The choice of the values of these 
parameters is dictated by phenomenology, and not determined by theory. 

Here we report a modified acceleration law which is essentially identical to that in STVG, but which 
comes from an entirely different line of reasoning. As motivation, we begin by stating an analogy 
with electrodynamics. If one considers a microscopic [corpuscular] medium of molecules in an external 
electromagnetic field, averaging of the Maxwell equations for the microscopic system results in the 
inclusion of polarization terms coming from the induction of moments on individual molecules. This is of 
course well known from electrodynamics. In an elegant analysis motivated by the study of propagation 
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of gravitational waves in a material medium, Szekeres [I] developed an analogous description for general 
relativity, via averaging of Bianchi identities. A modification results from quadrupole moments induced 
on 'molecules' which in the present context are galaxies, whose 'atoms' are stars. Zalaletdinov and 
collaborators [5] , [6] extended this analysis to the static weak-field Newtonian approximation, and showed 
that the effect of quadrupole polarization is to modify Poisson's equation to a fourth order biharmonic 
equation for the gravitation potential. In the present paper we obtain the modified acceleration law 
implied by this modification of Poisson's equation, and point out its significance for the observed galaxy 
rotation curves. 

The fact that our present analysis gives the same acceleration law as STVG is all the more intriguing 
because we were not looking for an explanation of STVG. Rather, it was a subsequent realization that 
the two laws are the same. Nonetheless, our presentation here emphasizes also the similarity between 
the two results - the commonality in spite of apparently different origins suggests the possibility that 
this may be a genuine effect, where the quadrupole induced modification can be mimicked by a field 
theoretic [fifth force] modification of general relativity. It should be emphasized that at this stage our 
analysis is intended to be illustrative, so as to bring out the physics involved, rather than to dwell into 
precise fitting of galaxy rotation curves. 

The plan of the paper is as follows. In Sections 2 and 3, we recall the acceleration law of STVG, 
and the derivation of the biharmonic equation for the gravitational potential. In Section 4 we obtain the 
acceleration law which results from this biharmonic equation and in Section 5 we compare it with the 
result from STVG. Section 6 briefly discusses the theoretical determination of the free parameters, ad 
Conclusions are given in Section 7. 



2 Galaxy Rotation Curves in the STVG Theory 



Here we briefly recall the form of the acceleration law, as derived in STVG theory. Details can be found 
in Q], and [3]. 

In the STVG theory, the action is given by 



S — Sgtscv + S<h + Ss + S_ 
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Here V(<p) is the potential for the massive vector field <f>^ having mass fi and coupling to; V(G),V(uj) 
and V(fJ,) are the potentials associated with the three scalar fields G(x), u)(x) and jJi{x), respectively, and 



dai>v - 9 v (j>u 



(5) 



By assuming that A = 0, V(<p) = 0, uj ad /i to be constants, and by considering the motion of a test 
particle coupled to gravity and to the vector field, it was shown that the law of acceleration is given by 



G X M exp(-r/r ) 
a(r) = J— + K 



r 

I + — 



where G M is defined to be the effective gravitational constant at infinity 

Goo = g(i - 




(6) 



(7) 
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and ro = 1/m- Here, Mo denotes a parameter that vanishes when u> — and G is Newton's gravitational 
constant. The constant K is assumed to equal 



K = G^MMq. 



(8) 



The choice of K, which determines the strength of the coupling of Bp V to matter and the magnitude of 
the Yukawa force modification of weak Newtonian gravity, is based on phenomenology and Mq is a free 
parameter of the theory. 

By using ([7]), one can rewrite the acceleration in the form 



a{r) 




1 - exp(-r/r ) 1 + 



0) 



It is assumed [and justified by further considerations] that one can generalize this to the case of a mass 
distribution by replacing the factor GM/r 2 in (J9j) by GM.(r)/r 2 . The rotational velocity of a star v c is 
obtained from v 2 (r)/r = a(r) and is given by 




1 - exp(-r/r ) 1 + 



ro 



1/2 



A good fit to a large number of galaxies has been achieved with the parameters: 
M = 9.60 x 10 11 M Q , r = 13.92 kpc = 4.30 x 10 22 cm. 



(10) 



(11) 



In the fitting of the galaxy rotation curves for both LSB and HSB galaxies, using photometric data to 
determine the mass distribution A4(r), only the mass-to-light ratio (M/L) is employed, once the values 
of Mo and ro are fixed universally for all LSB and HSB galaxies. Dwarf galaxies are also fitted with the 
parameters 

M = 2.40 x 10 u M Q , r = 6.96 kpc = 2.15 x 10 22 cm. (12) 

We will now see how a modified acceleration law of the form Eqn. ((9]) is obtained without modifying 
general relativity, as a result of considering polarization effects on the averaged gravitational field. As 
shown in the work of Szekeres and Zalaletdinov et al., which we recall in this next section, such a 
consideration leads to the Poisson equation being modified into a fourth order biharmonic equation. In 
the subsequent section, we derive the modified acceleration law by solving the biharmonic equation. 



3 Quadrupole gravitational polarization and the biharmonic equa- 
tion for the potential 

This section summarizes the work of Szekeres [J] and Zalaletdinov et al. [5], [5] leading to the derivation 
of the biharmonic equation. 

Thinking of galaxies as 'molecules' made of atoms [the stars], one would like to analyze how the 
averaged gravitational field inside a galaxy, is modified by the polarization of the molecules, due to 
the external pull of other galaxies. Indeed one has at the back of the mind the polarization induced 
modification of electromagnetic fields in a an electrically charged medium. 

In vacuum, if the Riemann tensor Rp, vpa is regarded as being the gravitational analog of the Maxwell 
tensor F^, the Bianchi identities 

R^u[pa-T\ = Rfivpc = (13) 

show striking resemblance to the "source-free" Maxwell equations 

iW] = °> F fr = °- ( 14 ) 

In the non-vacuum case, the Weyl tensor 

Cp, V p a = Rpvpa — 9p[pRa]v ~ Rp[p9a\v + T;R9^[p9u\a (15) 
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is the gravitational analog of the Maxwell tensor. The Bianchi identities now take the form, 

C'fivpa — &Jp,vpi (16) 

where 

J pup = J\nu] P = -(Tpfcv] - ^g P [f,T u] ). (17) 
Eqn. (fT5|) is comparable to the Maxwell equations 

Ffr = W„ (18) 

if C^vpa is interpreted as representing the free gravitational field and J^ vp is a kind of matter current. 
And the corresponding conservation equation is: 

J'U = 0- (19) 

The gravitational field created by a number of particles represented by a microscopic energy- momentum 
tensor ^( imcro ) } s given by Einstein's equations 



l/i a n^r - K M mlcro ') 
y •ep — ^Opg — — k 

where k = 8ttG/c 4 . The Einstein equations for the microscopic distribution of gravitational molecules 



9 ae r e0 -^ v r, v = -Kty— > (20) 



arc then : 



m^P-P-^b - |£(r ) - Sl (T a )]dr a (21) 
dn dr a dr a 

9 ae r^-^g^ v = -ntf molec) . (22) 

Averaging the left-hand side of the Einstein equations using Kauffman's method of molecular moments 
and by taking into account the expression [3] for the tensor of gravitational quadrupole polarization Q^ vpa 
in terms of the covariant gravitational quadrupole moment q% vpcr 

Q n Vpa = c -i(£ a j q^6 4 (x-y a )dT a ) (23) 
brings the averaged Einstein equations in the form: 

- \gp°}R {0) = -*{T«r e) + t£ w) + ( 24 ) 

where T^ W ^ is Isaacson's energy-momentum tensor of gravitational waves and Tll ree) is the energy- 
momentum tensor of molecules 

T^(x) = c-'E a j mo MM^-^( Ta )]dT . (25) 

A model of a static weak- field averaged gravitating medium with quadrupole gravitational polarization 
is developed next. The model is based on the following three assumptions. 

1. The energy-momentum tensor of the molecules T^ 00 ^ satisfies the standard conditions for the 
Newtonian approximation to hold : 

rp(free) ^(frcc) rp (free) ^(frce) ^(frce) _ ^(free) u v _ 2 ic,c\ 
J 00 ^ 1 ij ' J 00 1 0i ' J 00 ~ - L p,v U U — fit . [£V) 

2. The macroscopic metric tensor gjxJ is static, 



dt 



(o) 

= 0, (27) 
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3. The averaged macroscopic metric obeys the weak field approximation : tensor gjxJ , 

9^=V + e V. ( 28 ) 
In general relativity the conditions (|26l) - ((28|) for a microscopic energy-momentum tensor ^( mlcro ) anc [ 



for the metric tensor goo = ~ ^1 + <?oi — 0, g^ — Sij lead to the Newtonian limit of the Einstein 
equations (|2"U)) and to result in the Poisson equation for the Newtonian gravitational potential cf> = <fi(x a ) 

V 2 (j) = 47rG/i. (29) 

We are interested in the Newtonian limit of the macroscopic gravity equations (|24l ) under conditions 
(f2"6"| - (|28[) for the macroscopic tensor gfij 

9$ = -(l + Tl> 9^=1, 9^=1, $=1, 9$=0, ptu, (30) 



which leads to a generalization of the Poisson equation which incorporates the effect of gravitational 
quadrupole polarization. 

For the case of the static weak-field macroscopic medium Isaacson's energy-momentum tensor of 
gravitational waves vanishes 

T^ W) = 0. (31) 

The matter quadrupole polarization tensor is modelled by assuming, in analogy with the electrostatic 
case, that 

QiOjO — CgCiOjO (32) 

and the gravitational dielectric constant e g is modelled by [I] as 

1 mA 2 c 2 
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N. (33) 



Here, A is the typical linear dimension of a molecule, m is the typical mass of a molecule, ujq is a typical 
frequency of harmonically oscillating atoms in the molecule and N is the number density of the molecules 
in the medium. 
For the metric 



we get 



fl8? = -(l + 5). *2 } = 1, flS J =l, 9^=1, g$}=0 (34) 



R 00 = -4v 2 0, (35) 
&■ 

1 d 2 cj) 

Ru = -^-^—o, {i is n °t summed over here.), (36) 

R = gTRw (37) 

= g m R m +g u R^ (38) 



(— 1 + 2<j)/c 2 )( — ^ 2 4>) + ( — (* ls summed over here.) (39) 

2 



c 



2 V 2 </). (40) 



Terms only up to first order in </> have been considered. The last term in eauation (|24[) is Qfi pya - The 
only non-zero components of this tensor are: 

e d 2 6 

Q«yo = ' 9 R ^o = -^q-q- (41) 
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Substituting these values in the trace of Eqn. (IM|) , we get 



gTR^-^RgTg^ = -k^t^ + -c 2 g^Q^ ua ) (43) 

-R = - K (T+±c 2 g 00 Q^) (44) 

R = k(hc 2 + ^V 4 </>) (45) 

|lv 2 ^ = ^(^ + |VV) (46) 

V 2 = 4 7 rG Ai + ^^V 4 (47) 

y3 _ 27rGeg y4 = 4 



This is the biharmonic equation for the macroscopic gravitational potential. The origin of the fourth 
order term lies in the fact that the averaging procedure leads to the second derivative of the quadrupole 
polarization term in the averaged Einstein equations, and furthermore the polarization tensor is modelled 
as being proportional to the electric part of the Weyl tensor. 



4 Solution to the biharmonic equation for the macroscopic grav- 
itational potential 

As reviewed in the previous section, in the model [5], the equation for the macroscopic gravitational 
potential is the biharmonic equation 

V"V - k 2 V 2 (j) = ~<inGk 2 Li(r) (49) 

where k is a fundamental constant of the theory, having dimensions of inverse length, and defined by 
k 2 = c 2 /2irGe g . We are interested here in the modification of the radial dependence of the potential. 
Assuming a separation of variables, the radial part of this equation is given by : 

0"" + V" - k 2 cf>" - -fcV = -47rGfeVr) (50) 
r r 

where a prime denotes a derivative with respect to r. 

We find the series solution of this equation using the standard Frobenius method around the regular 
singular point r = 0. 

Case 1: The vacuum solution fi — [homogeneous equation] : 

In this case we get the following solution for the acceleration a = — V(f> 

1 + t) Wr + ( Co t) 6 ~2kV + ( Co + t) 2^ + ( Co t) - 1? (51) 



Since we know that k > 0, this implies that in Eqn. (|51[) . terms proportional to e kr — > oo as r — > oo 
which is unphysical. So we set the coefficients of e kr to zero. Thus 

Co + y \ = d = -kC . (52) 

Hence 

-kr -kr Q 

a = C -— + C 0TTT -^| (53) 
kr k z r z r z 



-$ + C o|J(l + M (54) 



G 



Eqn. (1541 is the solution for acceleration for vacuum [homogeneous] case. The similarity with the solution 
in the STVG case is immediately apparent. The constants Co and C2 can be related by the following 
reasoning : For kr -C 1 we assume Newton's law of gravitation to hold, so that C2 = GM+Co/k 2 = G^M 
where Goo = G[l + Co/k 2 MG}. For kr 3> 1 the exponential term can be ignored, and G m represents the 
effective gravitational constant at large distances. In the concluding section we will give an argument 
for estimating the physical value for the correction term proportional to Co- 

Case 2: /i(r) 7^ [Inhomogeneous Solution]: 

Since we know the solution for the homogeneous case (/i = 0), we can construct the solutions for the 
inhomogeneous case using the homogeneous solutions. As is well known, if y = au + bv is a solution of 



y" + P(x)y' + Q(x)y = 
where a and b are constants, we can find the solution of 

y" + P(x)y' + Q{x)y = R(x) 



in the form 
where 

where W is the Wronskian 
In our problem, 



y = A{x)u + B{x)v 
Mx) = - / — dx, B(x) = J —dx 



W — uv — uv. 



A(r) — 4nGk J r smh(kr)/i(r)dr + oi, B(r) — — inGk 2 J r cosh(kr) fidr + 0,2 
Hence inside a medium with a density profile ji{r) the acceleration is given by 



(55) 

(56) 
(57) 
(58) 
(59) 

(60) 



1 



AirGk cosh(fcr) 



AirGk sinh(fcr) f cosh(fcr) 
' r cosh(kr) /idr + ai 



r smh(kr)fidr 
sinh(fcr) 



a 2 - 



kr 



dr 1 



(61) 



The last two terms in the bracket in the above equation are same as in the homogeneous case. These 
two will be reduced to a\e~ kr jk 2 r 2 (\ + kr) for the same reasons as discussed before Eqn. (f54"|) . Hence 
Eqn. (H]) becomes 



1 



AnkG cosh(fcr) 



r smh(kr)/idr- 



AnGk sinh(fcr) 



r cosh(kr)fidr 



dr + 



-kr 



k 2 r 2 



.(1 + Jfer)-^ 



(62) 



Eqn. Q gives the acceleration for any given fi(r). One can easily read off special cases from here 
Case I : /i = : 



(63) 



which is same as equation (|54[) . as expected. 
Case II: fi =constant=/io : 



a = —0 = — - 



AttG^o 



r 2 k cosh(fcr) 



{ / r sinh(kr)dr} 



dr- 



AnGfiQ 



r 2 fcsinh(fcr) 



{ r cosh(kr)dr} 



-kr 



(64) 



or 
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47rG/Jo aie 
a = r 



k 2 



— kr n 



(65) 



Eqn. (|63|) gives acceleration for a medium of constant density. 



5 Comparison with the solution in Scalar-Tensor- Vector Grav- 
ity theory 

As noted earlier, in STVG the acceleration outside of a body is given by 



a{r) 



GooM . K 



-oexp(--) 1 + - 



'"o 



(66) 



This form matches exactly with that derived in (|54[) for the biharmonic equation. Comparing Eqn. ([5 
with Eqn. (|66p . 



C 2 = GocM = G 1 



Co 



M, 



& 



= K = G^MM~ 

k = i. 

ro 



(67) 

(68) 
(69) 



The free parameter ro of STVG theory is thus now determined from microscopic considerations which 
determine k from the gravitational dielectric constant. In the concluding section we will comment on 
the possible physical origin of K . 

Next, in STVG theory, the acceleration inside a spherical mass distribution takes the form [as can be 
shown with some justification] 



GM(r) J /Mb 



1 - exp(-r/r ) H 

ro 



(70) 



This is the form used to fit galaxy rotation curves. We write the above equation in the following form: 

a(r) = ^A(r) (71) 

where 



A(r) = 



GM 



1 



1 - exp(-r/r ) 1 + 



(72) 



Now one would like to check whether the biharmonic equation also gives the same form of acceleration 
inside a spherical mass distribution. In general this would be difficult to check; below we describe a way 
which suffices for our purpose. 

The solution of the biharmonic equation for any given mass distribution fj,(r) is given by Eqn. (j4]). 
Rewriting (HJ), 

i f r f 

AnkGr cosh(fcr) / r sinh(/cr)/zdr— 



4nGkr sinh(fcr) / rcosh(kr)/id 



dr + A(r) 



with 



ai = k 2 G\J MM , a 3 = G^MM + GM 



(73) 
(74) 
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we see that there is double integral over fj,{r) in above equation. We construct a function of the acceler- 
ation and its derivatives which is independent of integrals. Such a function is 



where 



X (r) = /3"(r) - fc 2 /3(r) 
P = ~[r 2 (a(r)-A(r))} 



From the solution of the biharmonic equation it can be shown that (as expected) 

X = A-nGk 2 r^,(r) 

We also find \(r) from STVG theory also using equation ([7T|) . Now 

asTVa(r) = a bh (r), 
=> XSTVc(r) = Xbh(r) 



(75) 
(76) 

(77) 

(78) 
(79) 



where "bh" stands for biharmonic equation. Equating xstvg and Xbh, we get an equation of the form 



P(r)M(r) + Q{r)M'(r) + R(r)M" (r) + S{r)M"'(r) + T(r) = 0. 



where 



P(r) 
Q(r) 

R(r) 

S(r) 

and T(r) 



G 



= G 



-r/r r; / 

M 



2 




2 2e' r/r ° 3e~ r/r ° e' r / r ° 2e~ r/r ° 



+ 



r 2 ro 



4tt 



rr 



G\ --- 

r r 

e -r/r oGM fM^ 



p—r/ro. M&. P ~ r l r a. Msl 
M e V M e V M 



'"() 



M(r) = J fir 2 dr, 



(80) 

(81) 
(82) 

(83) 

(84) 

(85) 

(86) 



Eqn. (|80p contains one integral over /i(r). Dividing eauation (|80[) by P(r) and then differentiating, we 
get a differential equation of the form, 



where 

U(r) 
V(r) 
W(r) 
and X(r) 



U{r)n(r) + V{r)p!{r) + W{r)n"(r) + X(r)//"(r) = 

(^e r ^"M^-7rr (r + r ) + M (r 2 - 6r r + 5r 2 )^) , 

2e r / ro M^^-r Q (2irr + r ) + M ( - 3r 2 + llr r + 2(-l + e r / r °)r 2 
e r/r °M^j-r Q (r + 3r ) + M (V 2 + (-4 + e r / r °)rr + 3(-l + e r / r °)r 2 



(87) 

(88) 
(89) 
(90) 
(91) 
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We do not get any term independent of fi(r) in the above equation. The solution of Eqn. ([57)1 will tell 
us for what value of n(r), STVG theory and the biharmonic equation give the same acceleration inside 
a mass distribution. But it is not necessary to solve (|57|) . Instead we take the observed density profile 
H(r) for a specific galaxy type and see whether it satisfies (l87l) . 
We assume the following density profile ji(r) 



fx(r) 



47rr 3 



/3M(r) 



where 



and 



M(r) = 4vr / dr'r' 2 n(r') 
Jo 

1 for HSB galaxies, 

2 for LSB & Dwarf galaxies. 



(92) 



(93) 



(94) 



We have taken the data set of ten LSB, HSB and dwarf galaxy masses from Table (3) in [2j- Then we 
plot the left hand side of Eqn. ([57]l 



STVG - BH 



M 2 



{U{r)fi(r) + V(r) f /{r) + W(r)n"{r) + A»,i'») 



(95) 



after making it dimensionless using parameters ro and Mq. We have made the assumption that k = 1/fo- 
We have taken data set of 10 LSB, HSB and dwarf galaxies and for each of them we have plotted 
STVG - BH vs. r for three ranges of r- (0.0001 - 1) kpc, (1 - 100) kpc and (100 - 500) kpc. The plots 
are shown in Fig. 1-10. We see that for the range 1 — 100 kpc STVG — BH is very small, of the order 
of 10~ 4 or 10~ 3 . And for other ranges STVG — BH takes high values. Since the range of length scales 
of interest in a galaxy lie between 1 — 100 kpc and STVG — BH is very small in this range, it can be 
said that the solution coming from STVG theory matches with the solution of the biharmonic equation, 
inside the medium, for the observed density profile. In this sense, the modified acceleration law resulting 
from the biharmonic equation can also explain the flatness of rotation curves, if k and Mq are assumed 
to take the values needed to explain observations. 





Figure 1: for galaxy of mass 0.13 x 10 Mq and r c — 0.53 kpc 



6 Determination of parameters from theory 

In the considerations based on the quadrupole polarization effect on the galactic field, k = 1/ro is a 
fundamental parameter derivable from microscopic considerations. We can find the value of k using the 
form given in [B] 

1 27I "G 

-rt (96) 



k 2 
where 



(97) 
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Figure 2: for galaxy of mass 0.42 x 1O 1O M0 and r c — 0.66 kpc. 
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Figure 3: for galaxy of mass 0.40 x 1O 1O M0 and r c — 0.82 kpc. 
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Figure 4: for galaxy of mass 0.38 x 1O 1O M0 and r c — 0.57 kpc. 
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Figure 5: for galaxy of mass 2.26 x 1O 1O M0 and r c — 1.06 kpc. 
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Figure 6: for galaxy of mass 3.08 x 1O 1O M0 and r c = 1.58 kpc. 






Figure 7: for galaxy of mass 5.46 x 10 10 M Q and r c = 1.40 kpc. 






Figure 8: for galaxy of mass 7.95 x 1O 1O M0 and r c — 1.36 kpc. 




Figure 9: for galaxy of mass 9.12 x 1O 1O M0 and r c = 1.04 kpc. 
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Figure 10: for galaxy of mass 20.19 x 1O 1O M and r c = 1.84 kpc. 



where wo=typical characteristic frequency of harmonically oscillating particles in the molecule, A — 
average linear dimension of a typical molecule, m = the average mass of a molecule, D is the mean 
distance between molecules, and N = average number of molecules per unit volume, k can also be 
written in the following form, 

' M^V- (98) 



1.2 



Here the dimensionless factor 9, 



reflects the nature of the field responsible for bounding discrete matter constituents into molecules. Here 
/i is the mass density of the medium. If 9 ss 1, the molecules of a self-gravitating macroscopic medium 
are considered to be gravitationally bound. 

In our calculation, the macroscopic medium is a galaxy cluster which is made up of galaxies ['molecules'] 
of which the stars are atoms. Since in galaxy clusters, galaxies are gravitationally bound, 
S«l, and assuming A w 100 kpc, and D«l Mpc. So we get 

q / in 5 \ 3 

1/fc2 " 8 X (iTT^) x( 100k P c ) 2 ( 10 °) 

or l/k s» 2 kpc . Thus it follows that that the estimated value of 1/k — ro is of the same order as desired 
in STVG theory. Of course this is only a ballpark estimate intended to show that k lies in the kpc range 
as opposed to being somewhere very far from the kpc range. 



7 Conclusions 

We found the general series solution to the biharmonic equation for the macroscopic gravitational po- 
tential. We saw that the form of acceleration is very similar to that in STVG pQ in vacuum. Then 
we investigated further if the form of acceleration matches inside the matter distribution also. What 
we find is that both the solutions match for the observed matter density profile of galaxies. Also, the 
value of the parameter ro in STVG [2] is of the order of k in the theory of macroscopic media in general 
relativity [BJ. 

While we have not offered an analogous estimate for Mo (equivalently K or Goo), it is clear that Mo 
also has to follow from microscopic considerations, since k is the only new parameter that appears in the 
biharmonic equation. A suggestive estimate may be obtained by rewriting the biharmonic equation by 
taking the fourth order term to the right : 

V2 ^ 4 ^ G ( 1+ 47gWm) ( 101 > 

The correction to G on the right hand is of the order 1 / k 2 r 2 which as we seen above is of order unity in the 
region of interest, thus suggesting that M is comparable to Mo, as assumed for explaining observations. 
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A few words about comparison with MOND. In MOND, the acceleration of a test particle is assumed 
to be given by the relation 

a = a N n( — ) (102) 

where is acceleration in Newtonian mechanics and is a new fundamental constant of nature. For 
small accelerations [large distances] it is assumed 

M— ) = — ■ (103) 

whereas fi approaches the value one for accelerations encountered in the solar neighborhood. In our work 
and in STVG theory as well, the acceleration is of the form 

a(r)=a N ^p-A(r) (104) 
which is also a form of modified Newtonian dynamics. But the extra factor 

M ( r ) A , x / a N , intn 
——A(r) ^ — . 105 
M ao 

The fact that the modified acceleration law discussed in this paper also arises from apparently indepen- 
dent considerations as in STVG, is intriguing and suggests that further investigation would be of interest. 
Especially, a better understanding of the physics involved, which goes beyond the mathematical treat- 
ment, is essential. Also, one would like to understand better the relevance and validity of the functional 
form for the gravitational dielectric constant used here. 

Of course no claim is made that we have proposed a definitive alternative to dark matter in galaxies, 
especially considering that any such proposals must also be studied for their role on cosmological scales, 
where dark matter is crucial for understanding structure formation. Our principal motive is to highlight 
that the role of averaged gravitational effects may not yet be fully understood on galactic scales, and 
even a percentage effect would be of interest and worth establishing accurately. 
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ADDENDUM: After the completion of this work two important points were brought to our attention 
by Viktor T. Toth 7 . Firstly, in a more recent treatment [5] than the one discussed by us above 
the parameters Mq and ro can, in fact, be obtained from three universal constants, which are fixed 
by observations ranging from the solar system to cosmological scales. It will be interesting to see in 
future work as to how the fundamental parameter k of the biharmonic equation relates to these universal 
constants. 

Secondly, in [5] it is shown in their Eqn. (65) that ro is proportional to \/M, where M is the source 
mass. Remarkably it turns out that in our case too this proportionality holds as is evident from our 
Eqns. and ([57| . (This is because the gravitational dielectric constant is proportional to the mass of 
the molecule) . This further strengthens the likelihood of a fundamental connection between STVG and 
our work. 

We would like to thank Viktor Toth for bringing these important observations to our attention. In 
further work we plan to explore the parallel between the two approaches in greater detail. 
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